Yukawa particles confined in a channel and subject to a periodic potential: ground 

state and normal modes 
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We consider a classical system of two-dimensional (2D) charged particles, which interact through 
a repulsive Yukawa potential exp{—r/X)/r, confined in a parabolic channel which limits the motion 
of the particles in the y-direction. Along the a;-direction, the particles are also subject to a periodic 
potential substrate. The ground state configurations and the normal mode spectra of the system 
are obtained as function of the periodicity and strength of the periodic potential (Vq), and density. 
An interesting set of tunable ground state configurations are found, with first and second order 
structural transitions between them. A magic configuration with particles aligned in each minimum 
of the periodic potential is obtained for Vo larger than some critical value which has a power law 
dependence on the density. The phonon spectrum of different configurations were also calculated. 
A localization of the modes into a small frequency interval is observed for a sufficient strength of 
the periodic potential. A tunable band-gap is found as a function of Vo- This model system can be 
viewed as a generalization of the Frenkel and Kontorova model. 

PACS numbers: 64.60.Cn, 82.70.Dd, 63.20.D- 



I. INTRODUCTION 

Two-dimensional system (2D) are often created in the 
presence of a substrate which may induce a periodic 
potential on the particles. In the pioneering experimental 
work of Chowdhury et al. Q, the authors studied a 2D 
colloidal system under influence of an one-dimensional 
(ID) periodic potential. An optical tweezer was used to 
trap the colloids by laser beams. For very high values of 
light intensity, crystallization of the colloidal suspension 
was observed, when the periodicity of the substrate (pe- 
riodic potential) was chosen to be commensurate to the 
mean particle distance. Laser induced freezing which is 
caused by the suppression of thermal fluctuations trans- 
verse to the ID periodic substrate was found (liquid-solid 
transition) jsj. The system studied in Ref. [H is related 
to the colloidal molecular crystal (CMC) and received 
attention recently due to important applications in pho- 
tonic and phononic crystals [1, Q . 

Specifically, CMC occurs when the number of colloids 
is an integer multiple of the number of substrate min- 
ima, and has been investigated using in simulations 0, Q 
and realized experimentally jTi]. CMC is an interesting 
experimental system to study order and dynamics in 2D 
since typical particle size and relaxation times permit, 
e.g, to use digital video-microscopy to track particle tra- 
jectories, allowing a deeper investigation of the physical 
behavior of the system [8|. 
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Originally, CMC was proposed for 2D system in 2D pe- 
riodic potential (substrate). As known the dimensional- 
ity of the system plays an important role in many physical 
properties of distinct physical phenomena. In this sense, 
an interesting question is how the ordered structures and 
physical properties would be influenced by the dimen- 
sionality of the periodic substrate. Recently, Herrera- 
Velarde and Priego d, [13] studied a 2D system of re- 
pulsive colloidal particles confined in a narrow channel 
and subject to an external ID periodic potential, which 
could be seen as the ID version of the CMC. The main 
focus of the study was the role of the substrate in the 
mechanisms that lead to a variety of commensurate and 
non-commensurate phases, its effect on the the single-file 
diffusion regime, and the pinning-depinning transition in 
ID systems. The ID character of the channel was rep- 
resented by a hard wall potential. Due to the repulsive 
interaction between the particles and the nature of the 
confinement potential the density across the channel was 
found to be non-uniform with a higher density at the 
channel edges. 

In the present paper we study the ordered configu- 
rations and the phonon spectrum of a 2D system of 
repulsive (Yukawa interaction) particles confined in a 
parabolic channel and subjected to a ID periodic po- 
tential alon g th e channel. As compared to the systems 
in Refs. [ol Il0l|. and due to the parabohc shape of the 
confinement an opposite density distribution is observed, 
with particles more concentrated at the central region of 
the channel. As shown previously for finite size clusters 
of repulsive particles, the confinement potential is deter- 
minant, e.g. for melting and phonon spectra [ll|. The 
ID parabolic confinement introduces a quasi-lT) (QID) 
character to the system in the sense that particles are still 
allowed to move freely in the perpendicular direction of 
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the confinement potential. 

The interplay between the repulsive inter-particle in- 
teraction and the periodic potential determines the dif- 
ferent ground state configurations. Our model system of 
Yukawa particles can be realized experimentally using: 
i) a dusty plasma ii) colloidal systems [l^ \Va^ 

and iii) electrons on liquid helium [13, HI]- A dusty 
plasma consists of interacting microscopic dust particles 
immersed in an electron-ion plasma. The dust particles 
acquire a net charge and the Coulomb interaction be- 
tween the dust-particles is shielded by the electron-ion 
plasma resulting in a Yukawa or screened Coulomb inter- 
particle interaction. The dust particles are confined to a 
two-dimensional layer through a combination of gravita- 
tional and electrical forces. By microstructuring a chan- 
nel in the bottom electrode of the discharge it is possible 
to laterally confine the dust particles as was realized in 
Refs. [UH^Ii the strength of the ID confinement po- 
tential can be varied by the width of the channel or the 
potential on the bottom electrode. When the width of 
the channel is microstructured into an oscillating func- 
tion along the channel it will result in a periodic potential 
along the channel. 

Alternatively, one can confine charged colloids, that 
move in a liquid environment containing countcrions, into 
microchannels as e.g. recently realized experimentally in 
[2^ 1 . In this case the inter-colloid interaction can be mod- 
eled by a screened Coulomb interaction and the confine- 
ment potential is a hard wall potential. By changing the 
depth profile of the micro-channel it has been shown in 
[25| that the confinement potential can be tuned into 
a harmonic potential. Micro-structuring the width of 
the channel into an oscillating function along the channel 
will result into an additional periodic potential along the 
channel. 

In a previous work [26| the ordered configurations of 
Yukawa particles confined to QID were studied. A phase 
diagram was obtained as function of the particle density 
and inverse Debye screening length which is a measure of 
the strength of the inter-particle interaction. The compe- 
tition between the lateral confinement and the screened 
Coulomb interaction resulted in different phases where 
the particles are ordered in chains. The most well studied 
phases are the one- and two-chain configurations where 
the transition between those two phases occurs through a 
zig-zag transition. The latter is a continuous transition as 
found theoretically for mono- [ 27.] and bi-disperse [28|, |29i| 
systems, and experiment ally |22| . |23| with a power law de- 
pendence on the width [23, [30| . Here we are interested 
to investigate how the phase diagram will be modified 
when an additional ID periodic potential is present. For 
example, how the zig-zag transition will be modified by 
the periodic potential. 

The present paper is organized as follows. In Sec. [Tl] 
we describe the model system and methods used in the 
calculation of the properties. In Sec. IIIII we present the 
results for the different ground state configurations. In 
Sec. IIVI the normal mode spectra for the one and two- 



chain regimes are presented for different intensities of the 
periodic potential. Our conclusions are given in Sec. [V] 



II. THE MODEL 

Our system consists of identical point-like particles in- 
teracting through a screened Coulomb potential. The 
particles are allowed to move in a two-dimensional (2D) 
plane and are subject to an external parabolic confine- 
ment in the y-direction and a periodic substrate potential 
along the x-direction. A sketch of the present model sys- 
tem is shown in Fig. [1] The total interaction energy of 
the system is given by 
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where e is the dielectric constant of the medium in which 
the particles are moving in, A is the Debye screening 
length, Vq is the strength of the periodic substrate poten- 
tial, L is the periodicity of the substrate potential, and 
r'i = is the position of the i*'' particle. In order 

to keep in Eq. ([l} only the parameters which rule the 
physics of the system, it is convenient to write the en- 
ergy and the distances in units of Eq — {mujQq'^ /2e^y^'^ 
and To — {2q^/m€UjQy^^, respectively, and the screening 
parameter k = tq/A, where. We also define the dimen- 
sionless strength of the substrate potential Vq = Vq/Eq 
and r = r' /vq. In so doing, the expression for the energy 
is reduced to 

H=E^-^+E^'^ + ^oE-s(^) (2) 
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As can be observed from Eq. ([2]) the system is a func- 
tion of the parameters, k, Vq, L, and the density. In our 
numerical calculations k = 1, which is a typical value 
for dusty plasma and colloidal systems. We introduce a 
distance ao, which is defined as the distance between par- 
ticles when Vb = 0. The density (n) is the ratio between 
the number of chains Nch and oq, i.e. n = Nch/o-o- In 
this case, the system self-organizes in a multichain-like 
structure pG^ . 

The model studied in this work is related to the 
Frenkel-Kontorova (FK) model, which is a simple one- 
dimensional model that describes the dynamics of a chain 
of particles interacting with nearest neighbors in the pres- 
ence of an external periodic potential. This model was 
initially introduced in the 1930s by Frenkel and Kon- 
torova and was subsequently reinvented independently by 
others, notably Frank and Van der Merwe [3l|. It pro- 
vides a simple and realistic description of commensurate- 
incommensurate transitions when thermal fluctuations 
are unimportant. 

In the present work the quasi-lD character of our sys- 
tem makes it different from the ID FK-model. The par- 
ticles have the additional freedom to move perpendicular 



3 



to the chain which leads to a rich set of new phases. 
The presence of two length scales in the FK-model, i.e. 
the inter-particle distance and the periodicity of the ID 
potential, is the reason of the complex behavior of the 
model. The inter-particle potential favors a uniform sep- 
aration between the particles, whereas the V(x) tends to 
pin the particles at the minima of the periodic poten- 
tial. This competition between both interactions is often 
called frustration or length scales competition. 




FIG. 1: (Color online) A sketch of the model system. 

The minimum energy configurations are obtained by 
numerical and analytical calculations. In the numerical 
simulations, we typically considered 100-200 particles, to- 
gether with periodic boundary conditions in the uncon- 
fined direction in order to mimic an infinite system. We 
do not consider friction in the present paper. In spite of 
the primary importance of friction to the motion of the 
particles in real systems, the ground state configurations 
are not aflFected by it. 

Notice that the substrate is defined in terms of the 
parameter L. Comparing L and ao, we define here an 
initially commensurate (IC) when {L/ao — p/q^ with p 
and q integers) and initially non-commensurate (INC) 
regime of the ordered structures when the ratio L/ao is 
a irrational number. It should be emphasized that in 
these cases the inter-particle distance oq is defined in the 
absence of a substrate (Vq = 0). In the case Vq ^ 0, 
it is expected that the mean distance between particles 
along a given chain a changes as a function of Vq , driving 
the system to new commensurate or non-commensurate 
configurations. 

III. GROUND STATE CONFIGURATIONS 

In this section, we present the results obtained ana- 
lytically and numerically for the ground state configura- 
tions. In the former, we calculate the energy per particle 
for different configurations as a function of the strength 
and periodicity of the substrate. We minimize such ex- 
pressions with respect to the different distances between 
particles. The configuration with lowest energy is the 
ground state. In order to predict which structures should 



be taken into account in the analytical approach, we also 
use molecular dynamic simulations as a complementary 
tool. The numerical method can give us some hints about 
which structures to consider. It should be noticed that 
one of the draw backs of the numerical technique is that 
in some cases there exist a larger number of meta-stable 
states, mainly in the limit of high densities where the 
system is found in a multi-chain structure. On the other 
hand, the numerical approach is the only way to obtain 
the ground state configurations in some incommensurate 
regimes, which will be analyzed in the next sections. 

We show here that depending on the periodicity of the 
substrate we can tune the ground state configuration, in- 
duce structural phase transitions and control the number 
of chains. This is interesting from an experimental point 
of view, since the number of chains can be associated 
with the porosity of the system, making it a controllable 
fiher. 

The main features of the present model system can be 
already seen in the more simple situations of the single- 
and two-chain regimes. For this reason, we limit our- 
selves to these cases, because it simplifies the physical 
interpretation of our results. 



A. Single-chain regime 

As an example, we study in this section systems with 
densities n — 0.5 and n = 0.7, which are found in the 
single-chain regime [1^ in the absence of the substrate 
Vb = 0. For n = 0.5 we consider the commensurate ratios 
L/flo = 1 and L/aQ = 2, while for n = \/2/2 we consider 
the non-commensurate regime with L/oq = V2. 

The simplest IC configuration is the trivial single-chain 
regime where each particle is sitted in a minimum of the 
periodic potential with n — 0.5 and i/op = 1. In this 
case, the configuration remains the same for any value of 
Vq. Notice that the cases in which L/ao = 1//, where 
/ > 1 is an integer, will exhibit the same behavior, since 
each particle is positioned exactly in a minimum of the 
substrate potential. On the other hand, the case L/a^ = 
I is very different and the particle configuration depends 
strongly on Vq, as will be shown in the next paragraph. 

In the IC case with n = 0.5 and L/oq = 2, for small 
values of Vb the particles are alternately located at the 
minimum and at the maximum of the substrate potential 
[see inset (a) in Fig. [2] . A sufficient increase of Vb forces 
the system to a new single-chain configuration in which 
a pair of particles is located at each minimum of the 
substrate [see inset (b) in Fig. [2]. A further increase 
of Vb pushes each pair of particles closer to each other, 
increasing the repulsive energy between them. 

For a critical value of Vb (~ 0.8) a structural transi- 
tion to the two-chain configuration is induced [see inset 
(c) of Fig. [2]. The system changes from a one to two- 
chains configuration. In the two-chain configuration the 
separation between particles in the x-direction of each 
minimum of the substrate is zero, which means particles 



are aligned along the y-direction. The separation dy be- 
tween chains does not change as function of Vq. In this 
particular configuration dy is ruled only by the compe- 
tition between the repulsive interaction between parti- 
cles and the parabolic confinement, being independent 
of the strength of the periodic potential. The type of 
transition observed here is different from the one found 
in Ref. [13] , where the authors demonstrated that in the 
absence of a periodic potential and in the presence of 
a parabolic confinement we have only continuous tran- 
sitions from one to two chains. In our case the one- to 
two-chains transition is clearly a first-order phase tran- 
sition as function of Vq. Notice that in the two-chain 
regime the system is re-organized in a final commensu- 
rate structure with a new ratio L/a — 1. We can define 
here a commensurate-commensurate transition between 
different orders of commensurability. 
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FIG. 2: (Color online) Nearest neighbor separation between 
particles in the x- (dx) and y-direction (dy) as a function 
of Vb for the case n — 0.5 and L/ao = 2. Three possible 
configurations are shown as insets. 

The expression for the energy per particle which is able 
to describe all phases observed in the case with n = 0.5 
and L/qq = 2 and in the case n = 1.0 and L/uq = 1 is 
given by 
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where Cx = dx/ L and Cy — dy/L are, respectively, the di- 
mensionless separations between particles within a min- 
imum of the periodic potential along the chain and per- 
pendicular to it. 
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FIG. 3: Ground state configurations for the case n = y/2 and 
L/ao = \/2 for: (a) Vo = 0.17 (b) Vo = 0.19. 



Now we discuss the INC regime with n = \/2 and 
L /flo = \/2- The same general behavior of previous cases 
can be observed here, with several structural transitions 
ruled by the strength of the periodic substrate Vq (Fig. 
[3]). For a large enough Vq the system can be found in a 
final commensurate regime with L/a ~ 1/2, but now in 
the three-chain configuration with particles almost uni- 
formly distributed over chains. 



B. Two-chain regime 

In this section we consider the system with n — 1.0, 
where a two-chain configuration is found as the ground 
state for Vb = 0. Differently from what was observed 
in the one-chain configuration {n — 0.5), when L/at) — 
1.0, the two chain configuration remains, but the internal 
structure depends on Vq. This is shown in Fig. |4l[a), 
where the relevant internal distances [Fig. Iljb)] for the 
arrangement are presented as function of Vq. 

For Vb = 0.16 the system changes from a staggered 
[dx ^ 0) to an aligned (dx — 0) two-chain configu- 
ration, through a second order (continuous) structural 
transition, characterized by a discontinuity in the second 
derivative of the energy with respect to Vq. Notice that 
in the case L/ap = 1.0, there are always two particles per 
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period of the substrate potential, and in such a commen- 
surate phase the system is always found in the two-chain 
regime. 
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FIG. 4: (Color online) (a) Inter-particle separation as func- 
tion of Vb for n = 1.0 and L/ao — 1.0. (b) A sketch of the 
two-chain configuration with the distances dx and dy indi- 
cated. 



Next, we consider the more interesting case with n = 
1.0 and L/ao — 2. When Vb is increased some unusual 
configurations appear as shown in Fig. [5] Initially the 
particles move in the x-direction towards the minima of 
the periodic potential and at the same time each chain 
starts to break up into two chains [Figs. [SJc)]. The 
transition found here is second-order. 

With further increase of Vq the two inner chains move 
towards each other [see Fig. (Hd)] and merge into a sin- 
gle chain in the center [see Figs. [5](d,e,f)]. The particles 
in the outer chains move towards the minimum of the 
periodic potential [see Figs. [5](d,e,f)]. With further in- 
crease of Vq the pair of particles in the middle chain are 
pushed closer to each other and finally form a row of four 
particles directed along the y-direction and positioned in 
the minimum of the periodic potential [see Figs. [5](g,h)]. 
The configurations presented in Fig. [S] indicates a tun- 
able porosity of the system as function of Vq. This is very 
convenient feature if the system is settled to be used as 
a filter or sieve, as pointed out in Refs. [s^ [111, where 
superparamagnetic colloidal particles were self-assembled 
in chain-like structures and used for separation of DNA 
molecules. 
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FIG. 5: Ground state configurations for the case n — 1.0 and 
L/ao = 2 for different values of Vq. 



The movement of the different particles in the x- and y- 
direction as function of Vq is summarized in Fig. [6l where 
structural transitions are indicated by vertical dashed 
lines. Three second order structural transitions are ob- 
served as function of Vq and the number of chains varies 
from 2-^4^-3-^-4-^4. 

Finally, we study the case with n = 1.5 and L/ag — 
1.5. Here the commensurate ratio changes according to 
Vq. Initially, for Vq = 0, the system is arranged in two 
chains [see Fig. EJa)], which are displaced with respect 
to each other over half the inter-particle distance in each 
chain. There are two particles per unit cell, which char- 
acterize an initially commensurate (IC) configuration. 

When Vo increases the system transits to a four-chain 
configuration through a second or first order structural 
transition, with the outer chains having twice as many 
particles as the inner chains [Fig. El^b)]. Alternatively, we 
can also view this configuration as two chains of triangles 
as indicated in the shadowed region in Fig. [7fb) . In this 
case d2 > and > d^ and the length of the unit cell 
is di = d2 + d^. There are six particles in the unit cell. 
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FIG. 6: (Color Online) (a) The lateral j/-position of chains 
for the case n = 1.0 and L/ao = 2. The vertical dotted lines 
represents the values of Vb where structural transitions occur, 
(b) The particle position in the i-direction as function of Vo- 



as in the case Vb = 0. 

With further increase of Vq, the y-distance dg between 
the internal chains goes to zero and the system changes to 
the three-chain configuration [Fig. [7{c)] with the same 
number of particles in each one, and the central chain 
shifted by a/2 along the x-direction with respect to the 
outer chains, which are aligned along the y-direction. No- 
tice that in this case there are only three particles per unit 
cell. This is interesting since the number of normal modes 
is now half the one observed for the configuration pre- 
sented in Fig. [71[b) , where the number of particles in the 
unit cell is six. The reduction of the allowed excitation 
modes is controlled by the strength of the periodic po- 
tential, and this can be used as an important feature for 
possible application in phononics. For Vb > 0.5, particles 
in different chains are all aligned along the y-direction 
and located in each minimum of the periodic substrate 
[Fig. El^d)]. The trajectories of the different particles in 
the channel as function of Vq is visualized in Fig. |S1 

Again, the relation between the periodicity of the sub- 
strate and the distance between particles is different from 
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FIG. 7: (Color online) Ground state configurations for dif- 
ferent values of Vb for n = 1.5 and L/ao ~ 1-5. In (b) the 
relevant distances for the analytical calculation of the energy 
are presented. 



the case Vq = 0, L/uq = 1. This is interesting since we 
can change the commensurability of the system by chang- 
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FIG. 8: (Color Online) Particle trajectories for different val- 
ues of Vo near the potential minimum for the case n = 1.5 
and L/ao = 1.5. 



ing only the strength of the substrate potential. 

As presented in Figs. (He), U and [5jh), for a critical 
value of Vq the present model system is found in the spe- 
cial configuration where the particles are aligned along 
the confinement direction. Such a y-aligned configura- 
tion (YAC) occurs if the condition L/ao = Pi where p 
is an integer (> 1), is satisfied. In this case, if N is the 
number of chains of the initial structure (Vq = 0), then 
we find that the number of particles aligned along the 
y-direction in each minimum of the substrate potential 
is N.p, which is also the number of chains. The critical 
value of Vq for which the YAC phase can be induced is 
obtained by adding the interaction between particles and 
confinement energies. A general expression for the YAC 
is given by: 
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FIG. 9: (Color online) (a) Critical value of Vb as function of 
density for the YAC phase, (b) The distance between the par- 
ticles along the y-direction as function of the density in the 
YAC phase. The inset shows a sketch of the general ground- 
state configuration with all relevant parameters. The red line 
in both figures is the linear fit. 



77 -kNpj/n 

V =— y 



Npy J 



+ 



2n ^ v-^ p-kNp/n^p+{p-qyc^ 



E EE' 

9=1 l=q+l j 



\J f + {p- q)'^c^ 



Np-l Np k(l-q)cNp/r, 



(5) 



E E 



A^V ^ ^ (l-q)c 

^ q=l l=q+l ^ ^' 

E'^ 



2c^Np 



1=1 

if N.p is an odd number. 

The critical value of Vq and the separation d between 
particles in each minimum are presented in Fig. 1^1 A 
sketch of the configuration in each minimum of the peri- 
odic substrate with all relevant parameters is also shown 
as insect in FiglSJb) . A power law dependence of Vq and 
d/L on the density is found. 
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IV. PHONON SPECTRUM 

Next, we analyze the Vo-dependence of the normal 
mode spectrum. We follow the standard harmonic ap- 
proximation and take into account the periodicity of the 
system in the unconfined direction (a;-axis). 

The present model system is a strictly 2D system where 
the number of particle in the unit cell and the number of 
degrees of freedom per unit cell determines the number 
of branches in the phonon spectrum. If I is the num- 
ber of particles per unit cell, there will be 21 branches in 
the phonon spectrum, from which half of those branches 
correspond to oscillations along the chain, i.e along the 
X axis we have longitudinal modes, while the others are 
associated with vibrations along the confinement direc- 
tion (y axis transverse modes). If the particles in the unit 
cell oscillate in-phase, the mode is dominantly acoustical, 
while the opposite out-of-phase oscillation corresponds to 
an optical mode. In general, a normal mode can be clas- 
sified in one of the following classes: longitudinal opti- 
cal (LO), longitudinal acoustical (LA), transverse optical 
(TO), or transverse acoustical (TA). 

In the harmonic approximation the normal modes are 
obtained by solving the system of equations 
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k 

FIG. 10: (Color online) The phonon spectrum for different 
values of Vo in the case n = 0.5 and L/ao = 1. 



The frequencies for the one-chain configuration when 
we have Vb ^ and two-chain configuration can be given 
by: 



D,. 



i)Qv.,i 



0, 



(6) 



= \l jiBi +B3± \lBl + AB^Bq - 2B1B3 + Bi + sub) 



where Qi/.i is the displacement of particle j from its equi- 
librium position in the v direction, ^ and v refer to the 

spatial coordinates x and y, 5^u,ij is the unit matrix and for the longitudinal modes, and by 
D^iy.ij is the dynamical matrix, defined by 



(10) 



(7) 



where u is an integer assigned to each unit cell. The force 
constants are given by 



, . ^ ^ ^ expj-K^y {x - x')^ + {y- y'Y) 

cl^f.A'^) = d^d, == , (8) 

^y[x - x'y + [y- y'Y 

with {x — x') — distance between particles along the x- 
axis and {y — y') = interchain distance with {x,y) and 
{x',y') the equilibrium positions of the particles in the 
unit cell, and 



0) 



(9) 



The phonon frequency is given in units of ujq / \/2. As 
an example, the complete dynamical matrix for the one- 
chain and the two-chain regime are given in Appendix. 

The frequencies for the one-chain configuration in the 
case Vq = 0, are given by uii = vCTj" for the acoustical 
branch and Wt = + A2 for the optical branch, Ai and 
A2 are defined in Appendix. 



Bo + Ba± 



ABrBr 



for the 
_Bi, i?2, -. 



2^2 B4 + Bi 
(11) 

transverse modes. The expressions for 
, and are given in Appendix. Here suh — 
8Vb7r^cos(7rCa;) is the term related to the periodic sub- 
strate. The wave number k for the one- and the two- 
chains regimes is in units of 2-k/L, where L is the length 
of the unit cell in the a;-direction. 

In Fig. Iior a). the phonon spectrum for the one-chain 
configuration is presented for different values of Vb, fixed 
density n = 0.5 and L/a^ = 1. In this case, there is 
one particle per unit cell located in each minimum of the 
substrate resulting only in one longitudinal mode and 
one transversal mode. The frequency of the longitudinal 
mode increases with increasing Vb, and there is a gap 
opening at A: = 0. The reason is that the periodic po- 
tential acts locally as a parabolic confinement potential 
V{x) ~ Vo^x^ with frequency u = ^V^jm'^. The 
A: = gap corresponds with this frequency for not too 
small values of Vq. The transversal mode corresponds 
to particle oscillations in the y-direction and is therefore 
practically independent of Vb- 

In Fig. 111! the dispersion curves for n = 0.5 and 
L/oo = 2 are presented for different values of Vb. As ob- 
served in Fig. [21 the presence of the substrate (Vb ^ 0) 
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FIG. 11: The phonon spectrum for different values of Vo in 
the case n — 0.5 and L/ao = 2. 



modifies the number of particles in the unit cell in or- 
der that the number of branches of the phonon spec- 
trum is increased as compared to the case Vo = [Fig. 
niT a)]. For Vb 7^ there are two particles per unit cell 
and consequently four branches in the phonon spectrum. 
As Vq increases, the frequency of the LA mode also in- 
creases, which can be explained keeping in mind that for 
low values of Vq, there is a small electrostatic repulsion 
between neighboring particles, in order that particles os- 
cillate horizontally without major difficulties. The oppo- 
site behavior is found for the TO mode, i.e., decreases 
with increasing Vo. The distance between adjacent par- 
ticles in the same substrate minimum becomes smaller, 
and the repulsive force between them increases and acts 
as a retarding force. 

The LO mode has a rather different behavior as com- 
pared to the TO mode, i.e., there is a hardening of its 
frequency when Vq increases, which is a consequence of 
the larger repulsion due to the closer proximity between 
particles. For a sufficiently strong Vq [ Fig. [TlTc)] the 
normal mode spectrum becomes discrete, i.e. frequen- 
cies become independent on fc, which means the group 
velocity is zero and the modes become localized. 

As commented previously, for Vq > 0.8 there is a struc- 
tural phase transition to the two-chain configuration with 
particles aligned along the y-direction (YAC phase) in 
each minima of the substrate [Fig. [2] . Again, due to the 
strong confinement the modes are localized, and a dis- 
crete spectrum is found [Fig. [TlTc')] with large frequen- 
cies for the longitudinal modes. The transversal modes 
also present an increase in the frequency due the large 
confinement in the y-direction. 

Now we discuss the dispersion curves for the system 
with density n — 1.0 and L/gq = 1, [Fig. [12]. As pre- 
sented earlier [Fig. 0] , particles remain in the two-chains 
configuration for all Vq with changes only in the internal 
structure. Again the substrate potential induces gaps in 
the normal mode frequencies as presented in Fig. 1121 
The TA, TO, LA and LO modes increase with increasing 
Vo. 

In the case of the LA mode, for low values of Vq par- 
ticles are not aligned, having more freedom to oscillate 
in the horizontal direction. When Vq increases, the elec- 
trostatic force becomes larger (particles are now aligned) 
making oscillations along the channel harder. 

The LO mode also increases with increasing Vq- This is 
a consequence of the strength of the substrate potential, 
which trap particles in their equilibrium positions, reduc- 
ing the out phase oscillations of the particles. The TO 
frequency increases slightly, since the out-of-phase mo- 
tion is more difficult to occur. The TA frequency branch 
is almost independent of Vq because it corresponds to os- 
cillations in the y-direction and is therefore determined 
by the harmonic confinement potential with frequency 
Wo- 

For the YAC phase, Vq > 0.16, the normal mode 
spectrum becomes discrete [Fig. [T^ d)]. as in the case 
n = 0.5, L/a = 2.0. The modes are almost constant 
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FIG. 12: The phonon spectrum for different values of Vb in 
the case n — 1 and L/ao = 1. 



due to the strong confinement potential imposed by the 
substrate and the harmonic trap. 



V. CONCLUSIONS 

We investigated the structural and dynamical prop- 
erties of a two-dimensional system of repulsive particles 
confined by a parabolic channel and submitted to a one- 
dimensional periodic potential (substrate). The ground 
state configurations were obtained analytically and nu- 
merically, where for the latter we used molecular dynam- 
ics simulations. The phonon spectrum were also calcu- 
lated analytically for the one- and two-chain configura- 
tions through the harmonic approximation. 

The main features of the structure and normal mode 
spectrum were studied (for different densities) as a func- 
tion of the periodicity (L) and strength (Vq) of the sub- 
strate, which are experimentally tunable parameters in 
systems like e.g. colloids in the presence of a periodic 
light field composed of two interfering laser beams. An 
interesting set of ground state configurations with con- 
trollable porosity is observed mainly as a function of Vq, 
through several first or second order structural transi- 
tions. The structures are mainly ruled by the fact that 
particles tend to go to the minima of the periodic sub- 
strate, modifying the symmetry of the ordered structures. 
However, for small Vq the inter particle repulsive interac- 
tion dominates and the particles are found over all possi- 
ble positions in the periodic potential, including regions 
near to the maxima. For large Vq, particles are more and 
more attracted to the wells of the periodic potential. 

For some specific cases we found structural transitions 
where the number of particles in the unit cell of the peri- 
odic system is changed, implying e.g. a different number 
of branches in the phonon spectrum, which is an inter- 
esting aspect of the dynamical behavior of the system, 
specially for applications in phononics. 

The normal mode frequencies depend on the linear 
density of the system, periodicity and strength of the pe- 
riodic substrate. We observed gaps in the phonon spec- 
trum, which indicate that there are frequencies blocked 
by the crystal. For Vq beyond a critical value and for spe- 
cific values of the ratio L/ao the system is found in a spe- 
cial configuration were particles are aligned in each min- 
imum of the periodic substrate and perpendicular to the 
x-direction. For such a configuration the normal mode 
frequencies become independent of the wave vector and 
the modes localize into a small frequency of interval. 
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Appendix 

The matrix w^I — D (where I is the unit matrix and D 

is the dynamical matrix) is used in the calculation of the 
normal modes for the one- and two-chains configurations. 
The dynamical matrix for one chain configuration when 
Vo=0 is: 

io'^-Ai 

{0J^-LJ^)-A2 

where the quantities Ai and A2 are given by: 



i=i 

00 



-2Kri/n 



(2ri)3 



3cf, 4K^cf, 6Kcf, 



nri 



2k j 
n 



[I — exp{ikjL)] 



(A.6) 



'^J/" 2Kj K^f 



J- 



+ 



i=i 

+ VoTrn^cos{Trj) 



A2 = J2n^ 



-Kj/n 



1 + 



KJ 



[1 — cos{knj)] 

(A.l) 

[1 - cosikirj)] (A.2) 



The dimensionless wave number k is in units of 27r/Z/. 
The dynamical matrix to one-chain Vb ^ and two- 
chains configuration is: 



^2 -Bi- sub -B5 

Alu^-B2 -Be 

-B7 uj^ -B3- sub 

-Bg Aw^-Bi 



where Aw^ = — Wg- The quantities Bi, B2, -B3, -B4, 
and Bq are given by: 



B5 = J2^' 



°° g-2Kr/7^ 



(2r)3 



0-Cx ' H + — + — 



(A.7) 



(1 H ) [exp{ikL{j - c^))] 

n 



5e = E^ 



-Ikt jn 



(2r 



1 3 



3c?, 



X [exp{ikL{j - c^))], 



nr n . 



(A.8) 



3=1 



(2r)3 



(i-c.)'( 



3 6k 4k2n 
nr n"' ^ 



-(! + — )J+E"-^- 

X [1 — exp{ikjL)] 



-2Kj /n 



(2j) 



;V5 



4Kj , (2Kj)2 



2+ — + 



(A.3) 



5r = E' 



g— 2Kri/n 



(2r)3 



2K;ri 

(1 + ) [exp{ikL{j + c^))] 

n 



{j-l + c^f [- + — + — ) 



(A.9) 



B2 = Y.'. 



00 

E- 



^—1nT In 

?- 

(2r)3 
(2jr 



3c2 



1 + 



4k"^c; 



2„2 



6kc,, 



nr 



2Kj 



[1 - exp{ikjL)], 



(A.4) 



58 = E^ 



-2K.ri / n 



(2ri)3 



X [ea;p(iA;iy(j 



3c2 



4k"' c: 



2^2 



nri n 



(A.10) 



^3 = E"^ 



-(1 + 



g-2Kri/n 

(2r)3 



/ ■ 1 \2/' ^ 6k 4k2\ 
vr^ nr n^ / 



2Kri 



■DC' 

)]+E^ 



-2K,j /n 



2 + 



4kj ^ (2«:j)2 



n 



X [1 — exp{ikjL)] 



(A.5) 



where r = ^ {i - c^)^ + c^, n = ^/(i- 1 -h c^)^ -h c^, 

the dimensionless wave number k is in units of 2tt/L, 
i = V— 1 and sm& = 8Vb7r2cos(7rCx). 
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